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Abstract 

We investigate spatially flat isotropic cosmological models which contain a scalar field with an ex- 
ponential potential and a perfect fluid with a linear equation of state. We include an interaction term, 
through which the energy of the scalar field is transferred to the matter fields, consistent with a term that 
arises from scalar-tensor theory under a conformal transformation and field redefinition. The governing 
ordinary differential equations reduce to a dynamical system when appropriate normalized variables are 
defined. We analyse the dynamical system and find that the interaction term can significantly affect the 
qualitative behaviour of the models. The late-time behaviour of these models may be of cosmological in- 
terest. In particular, for a specific range of values for the model parameters there are late-time attracting 
solutions, corresponding to a novel attracting equilibrium point, which are inflationary and in which the 
scalar field's energy-density remains a fixed fraction of the matter field's energy density. These scalar field 
models may be of interest as late-time cosmologies, particularly in view of the recent observations of the 
current accelerated cosmic expansion. For appropriate values of the interaction coupling parameter, this 
equilibrium point is an attracting focus, and hence as inflating solutions approach this late-time attractor 
the scalar field oscillates. Ifence these models may also be of importance in the study of inflation in the 
early universe. 
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1 Introduction 



A variety of theories of fundamental physics predict the existence of scalar fields , motivating the study 
of the dynamical properties of scalar fields in cosmology. Indeed, scalar field cosmological models are of great 
importance in the study of the early universe, particularly in the investigation of inflation [^Q. Recently 
there has also been great interest in the late-time evolution of scalar field models. 'Quintessential' scalar 
field models (or slowly decaying cosmological constant models) |^,|^ give rise to a residual scalar field which 
contributes to the present energy-density of the universe that may alleviate the dark matter problem and can 
predict an effective cosmological constant which is consistent with observations of the present accelerated 
cosmic expansion 0,1). 

Models with a self-interaction potential with an exponential dependence on the scalar field of the form 

V = Ae'^'^, (1) 

where A and k are positive constants, have been the subject of much interest and arise naturally from theories 
of gravity such as scalar-tensor theories or string theories |3| . Recently, it has been argued that a scalar field 
with an exponential potential is a strong candidate for dark matter in spiral galaxies |^ and is consistent 
with observations of current accelerated expansion of the universe . 

A number of authors have studied scalar field cosmological models with an exponential potential within 
general relativity. Homogeneous and isotropic Friedmann-Robertson- Walker (FRW) models were studied by 
Halliwell ||ll[] using phase-plane methods. Homogeneous but anisotropic models of Bianchi types I and HI 
(and Kantowski-Sachs models) were studied by Burd and Barrow Bianchi type I models were studied by 
Lidsey ||l^ and Aguirregabiria et al. , and Bianchi models of types HI and VI were studied by Feinstein 
and Ibahez JTsf . A qualitative analysis of Bianchi models with fc^ < 2 (including standard matter satisfying 
standard energy conditions) was completed by Kitada and Maeda . The governing differential equations 
in spatially homogeneous Bianchi cosmologies containing a scalar field with an exponential potential reduce 
to a dynamical system when appropriate expansion- normalized variables are defined. This dynamical system 
was studied in detail in [ pT| (where matter terms were not considered). 

One particular solution that is of great interest is the flat, isotropic power-law inflationary solution which 
occurs for fc^ < 2. This power-law inflationary solution is known to be an attractor for all initially expanding 
Bianchi models (except a subclass of the Bianchi type IX models which will recoUapse) [^|l^. Therefore, 
all of these models inflate forever; there is no exit from inflation and no conditions for conventional reheating. 

Recently cosmological models which contain both a scalar field with an exponential potential and a 
barotropic perfect fluid with an equation of state 

P = (7 - l)Ai, (2) 

where 7 is in the physically relevant range 2/3 < 7 < 2, have come under heavy analysis. One class of 
exact solutions found for these models has the property that the energy density due to the scalar field is 
proportional to the energy density of the perfect fluid, and hence these models have been labelled matter 
scaling cosmologies [p"8|-p0| . These matter scaling solutions are spatially flat isotropic models and are known 
to be late-time attractors (i.e., stable) in the subclass of flat isotropic models |^"8|-p0| and are clearly of 
physical interest. In addition to the matter scaling solutions, curvature scaling solutions [pT| and anisotropic 
scaling solutions |^ are also possible. A comprehensive analysis of spatially homogeneous models with a 
perfect fluid and a scalar fleld with an exponential potential has recently been undertaken . 

Within the standard model of inflation (with, for example, a quadratic or quartic self-interaction poten- 
tial), in the early universe the microphysics of the scalar field leads to an accelerated expansion essentially 
driven by the potential energy (or vacuum energy) that arises when the scalar field is displaced from its po- 
tential energy minimum. Provided that the potential is sufficiently flat, the universe undergoes many e- folds 
of expansion and the matter content is driven to zero 1 24 - 28) . As the scalar field nears its minimum, the 



vacuum energy is converted to coherent oscillations of the scalar field (which corresponds to non-relativistic 
scalar particles). Eventually, since the scalar field is coupled to other (fermionic and bosonic) matter fields, 
these particles decay into lighter particles and their thermalization results in the 'reheating' of the universe 
(thereby accounting for the current entropy in the universe today). 
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Although the exponential models are interesting models for a variety of reasons, they have some shortcom- 
ings as inflationary models. While Bianchi models generically asymptote towards the power-law inflationary 
model in which the matter terms are driven to zero for k'^ < 2, there is no graceful exit from this inflationary 
phase. Furthermore, the scalar field cannot oscillate and so reheating cannot occur by the conventional 
scenario. Clearly, these models need to be augmented in an attempt to alleviate these problems. For exam- 
ple, exponential potentials are only believed to be an approximation, and so the theory could include more 
complicated potentials (although it is not clear what these other potentials should be). 

The goal of this paper is to examine how interaction terms, through which the energy of the scalar field 
is transferred to the matter fields, can affect the qualitative behaviour of models containing matter and a 
scalar field with an exponential potential. Such an interaction term, denoted by 5, arises in the conservation 
equations, via 

^{<j) + 3F0 + kV^ = -5 (3a) 
/i + i-iHii = 5. (3b) 

The form of such a term has been discussed in the literature within the context of inflation and reheating. 
An alternative to the conventional reheating model, in which the scalar field's energy is transferred to the 



matter due to scalar field oscillations, is the warm inflationary model |20,^,|32l in which an interaction term 
is significant throughout the inflationary regime (not just after slow-roll) and so the energy of the scalar field 
is continually transferred to the matter content throughout inflation and the matter content is not driven to 
zero. 

Several examples of interaction terms appear in the literature for models with a variety of self-interaction 
potentials. In particular, potentials which have a global minimum have attracted much attention. For 
example, Albrecht et al. |^l[ considered 5 = a(j)'^4>^~'^'^ (where a is a constant), derived from dimensional 
arguments, in the reheating context after inflation (with potentials derived from Georgi-Glashow SU(5) 
models). In a similar context, Berera considered interaction terms of the form 5 = aij)^ . Quadratic 
potentials and interaction terms of the form 5 = a(p^ and 5 = 00^0^ were considered by de Oliveira and 
Ramos and a graceful exit from inflation was demonstrated numerically. Similarly, Yokoyama et al. | |3^ 
showed that an interaction term, which is negligible during the slow-roll inflationary phase, dominates at 
the end of inflation when the scalar field is oscillating about its minimum; during this reheating phase it is 
assumed that the energy transferred from the scalar field is solely converted into particles. 

Within the context of exponential potentials, Yokoyama and Maeda and Wands et al. ||l^ considered 
interactions of the form 6 = a\/V. The main goal in both these papers was to show that power-law inflation 
can occur for k^ > 2, thereby showing that inflation can exist for steeper potentials. The main motivation for 
this work is the fact that exponential potentials which arise naturally from other theories, such as supergravity 
or superstring models, typically have k^ > 2. Wetterich |36[ considered interaction terms containing a matter 
dependence, namely 5 — a(j)ii, in which perturbation analysis showed that the matter scaling solutions were 
stable solutions when such interaction terms are included. In | ]36[ |, it was shown that the age of the universe 
is older when S is included and that the scalar field can still significantly contribute to the energy density 
of models at late times. In [||, certain string theories in which the energy sources are separately conserved 
in the Jordan frame naturally lead to interaction terms in the Einstein frame, although this is not specific 
to string cosmologies; any scalar-tensor theory with matter terms and a power-law potential will yield the 
same results |3^ . 

The scalar field in scalar-tensor theory can be related to the scalar field in the general relativistic scalar 
field theory via a conformal transformation. If the matter terms are separately conserved in the scalar-tensor 
theory (Jordan frame) then they lead to the following conservation equation in the Einstein Frame 

3g0 + kv) ^ ~ 0^ (4) 

' 2 

(and a similar equation for /i), where uj is the coupling parameter in the scalar -tensor theory. This therefore 
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leads to an interaction term of the form S = 



—acjjii, where 



(4-37) ... 
a= ^ (5) 

^ 2 

Similarly, in string theory {uj = —1), a similar interaction term arises in the Einstein frame which depends 
on the energy density of the matter (axion) field and a is of order unity. Finally, a term of the form 6 = afiH 
might be motivated by analogy with dissipation. For example, a fluid with bulk viscosity may give rise to a 
term of this form in the conservation equation [ p7| . 

If S does not depend on the matter energy-density, then unphysical situations may arise. For conven- 
tional interaction terms without ^, we have found numerically that /i is driven to zero in a finite time and 
subsequently the matter energy-density becomes negative (see for details). Thus, we shall include a factor 
fi in the interaction term in order to ensure that /i > 0. In it was shown that if (5 < at the equilibrium 
points then static solutions (i.e. (/) = 0) are not possible. Furthermore, we require that the sign of S be 
positive at equilibrium points representing inflationary phases, otherwise the matter fields will be "feeding" 
the scalar field and will redshift to zero even faster than in the absence of the interaction terms. 

This paper examines interaction terms of the general form 5 — SfiH (where S = 5{(p, V, H)) in the context 
of flat FRW models with emphasis in determining the asymptotic properties of these models. In particular, 
it will be determined whether these models can asymptote towards inflationary models in which the matter 
terms are not driven to zero. This would partially alleviate the need for reheating; since the matter content 
tracks that of scalar fleld (in this context) it is never driven to zero (unless both are driven to zero in 
which case the solution is not inflationary). However, a more comprehensive reheating model would still 
be necessary. The structure of the paper is as follows. In section ^, the governing equations are defined 
and the case 6 — studied in |Q is reviewed. In section ||, the case S = acfifi is studied, motivated by the 
conformal relationships between the Jordan and Einstein frame in string theory, and extends the work of 
p6t . In section ^ we shall discuss the results of the analysis and in section || we shall examine an interaction 
term of the form 6 = ajiH. The paper ends with conclusions in section |6[ 

2 Governing Equations 

The governing field equations are given by the equations (^) and 

iir=-i(7M + 0'), (6) 

subject to the Friedmann constraint 

H'^lif^ + l^' + V), (7) 

(an overdot denotes ordinary differentiation with respect to time t). Note that the total energy density of 
the scalar field is given by /i^ = ^ijp + V. The deceleration parameter for this system is given by 
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and is independent of the interaction term. 
Defining 



and the new logarithmic time variable t by 



V6H ' V^H 



dT 



(8) 



(9) 
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the governing differential equations can be written as the plane-autonomous system: 



2 2\ 

X -y ) 



2^ 



^kx + + 7(1 ~x^ - y^) 



(11a) 
(lib) 



where a prime denotes differentiation with respect to t ||^. Note that y = is an invariant set, corresponding 
to y = 0. The equations are invariant under y —y and t ^ —t and so the region y < is a time-reversed 
mirror of the region y > 0; therefore, only y > will be considered. Similarly, only fc > will be considered 
since the equations for the interaction terms to be considered are invariant under k — > —k and x — s- ~x. 
Equation (^ can be written as 



where 



n^^J^ = x^ 



V 



(12) 



which implies that < < 1 for > so that the phase-space is bounded. The deceleration parameter 

is now written 

9 = -1 + 3x2 + ^7 (1-^2 -2/2). (13) 



2.1 Comments on arbitrary 5 



The fact that q is independent of the interaction term implies that the region of phase space which represents 
inflationary models is the same for all of the models considered. Namely, g = occurs along the ellipse 
72/2 _ (2 — 7)a;2 -|- i(37 — 2). For any value of 7, the lines intersect the boundary of the phase space at 
a;2 = i 

It is possible to make some qualitative comments about the system ( pi] ) for arbitrary 5. First, the location 
of equilibrium points on the boundary x^ + y^ = 1 (f2 = 0) are independent of the choice of such interaction 
terms; the three equilibrium points (and their associated eigenvalues) which exist on the boundary for any 
S are: 



IC+ : 


{x,y) = (4 


-1,0) 




(Ai,A2) = 


(^3(2 


IC- : 


ix,y) = {- 


-1,0) 




(Ai,A2) = 


(^3(2 


Ps ■■ 




k 




(Ai,A2) = 


H 



3 1^=+' V 2 



3 I'c- 



[6 - k^] 



37 — k~ 



V6 + k 



V6- fc 



3 



(14a) 



(14b) 



(14c) 



The points K.^ represent the isotropic subcases of Jacobs' Bianchi type I solutions |3^ (subcases of the 
Kasner models), generalized to include a massless scalar field. These solutions are non-inflationary {q = 2). 
The point Pg, which exists only for fc2 < 6, represents the FRW power-law model [|ri|,|T^ and is inflationary 
for fc2 < 2 (g = i [k^ " 2] ) . Although these three points exists for any S, the interaction term does affect the 
stability of these solutions, as is evident from the eigenvalues in (|l4|). In particular, for k^ < 37 the point 
Ps can become a saddle point if 

>3(37"fc2). 



5 



Hence, if the interaction term is significant, solutions will spend an indefinite period of time near this power- 
law inflationary model, but will then evolve away and typically be attracted to another equilibrium point in 
some other region of the phase space. 

Matter scaling solutions (i.e. those solutions in which 7^ = 7) , denoted by Ts in , exist only in special 
circumstances when such interaction terms are present, and occur at the point 



2 k' V 2fc2 



^y's = - V v^s = V — 7rr3 — ■ (15) 



6\^. (16a) 



Substituting these solutions into equations (|ll|) yields 

^, ^ (37 -fc^) 

3\/6fc7 "'^5' 

y' = 0. (16b) 

Hence, the matter scaling solutions will be represented by an equilibrium point only if ^|^^ — (or in the 
special case 87 = fc^, which is typically a bifurcation value). For the simple forms for d given in the literature 
and those used in this paper, this condition will not be satisfied and so the matter scaling solutions cannot 
be asymptotic attracting solutions. 

However, an analogous situation does arise. In particular, any equilibrium point within the boundary of 
the phase space will satisfy 

yl - ^^^4, (17) 

in which the scalar field is equivalent to a perfect fluid of the form — (7^ — l)/i^, but where 7^ 7^ 7. 
Consequently, any attracting equilibrium point within the phase space will represent models in which neither 
the matter field nor the scalar field is negligible and the scalar field mimics a barotropic fluid different from 
the matter field and therefore could still constitute a possible dark matter candidate. 

Finally, it can be shown that any equilibrium point within (but not on) the boundary will occur for 



xq < 0. For 2/ 7^ and Q ^ equation ( lib ), which does not depend on (5, yields 



7y' = 7+y^/ca;+(2-7)x^ (18) 

a relationship any such equilibrium point must satisfy. Now, 7(1/^ + x'^) < 7 since + x^ < 1 inside the 
boundary, and hence (|l8|) yields 

x(^x + ^^<0, (19) 
which cannot be satisfied for a; > (since A: > 0). 



2.2 Review of the case 5 = 

Copeland et al. ||3^ performed a phase-plane analysis of the system (^ij) for (5 = 0, and found five equilibrium 
points. One of the equilibrium points (denoted here by P) represents a flat, non-inflating FRW model [Q, 
for which = 1. For 2/3 < 7 < 2 this point is a saddle in the phase space. The flat FRW matter scaling 
solution {Fs) was found to exist for fc^ > 87 and was shown to be a sink. The equilibrium point /C"*" was 
shown to be a source for all k and /C~ a source for fc^ < 6. The FRW power-law model (P5) was shown to 
be a sink for fc^ < 87, and was shown to represent an inflationary model for < 2. The results found in 
psf are summarized in table [|. 



3 Interaction term of the form S = —a(f)ii 

In this section, an interaction term of the form 6 = ~a(f)fi (and hence 6 = —a<j)/H) shall be considered. 
Again, it will be assumed that fc > 0. The explicit sign choice for 5, with the assumption that a > 0, is to 
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< fc^ < 2 


2<k 


^ < 6 


A:^ > 6 






< 37 


k'^ > 37 




p 


saddle 


saddle 


saddle 




(NI) 


(NI) 


(NI) 




source 


source 


source 




(NI) 


(NI) 


(NI) 


K.- 


source 


source 


saddle 




(NI) 


(NI) 


(NI) 


Ps 


sink 


sink 


saddle 


DNE 




(I) (for fc2 < 2) 


(NI) 


(NI) 






DNE 


DNE 


sink 


sink 








(NI) 


(NI) 



Table 1: The equilibrium points for 6 = and their stability for various values of k. The label "(NI)" 
denotes non-inflationary models whereas "(I)" represents inflationary models, and "DNE" is used when an 
equilibrium point does not exist. 



guarantee that all equilibrium points within the phase space will represent models in which energy is being 
transferred from the scalar field to the perfect fluid, since it was shown that all equilibrium points within the 
phase space occur for a; < (0 < 0). Indeed, this is even true for the equilibrium point on the boundary 
= 0, since it is located at a; < 0. With this particular choice for S, equations (O) become 



-3a;(l - x^) 



y 



2^ 



^kx + + 7(1 - x^ ~ y^) 



There are five equilibrium points for this system: 

1. /C+: (a:,2/) = (+1,0), 17 = 0, q = 2. 
The eigenvalues for this equilibrium point are 

(Ai,A2) = (3(2-7)- V6a, 



3 r 



\/6 + /c 



(20a) 
(20b) 



(21) 



This equilibrium point is a source for a < y ^^(2 — 7) and a saddle otherwise 

2.K-: (a;,y) = (-1,0), = 0, q = 2. 

The eigenvalues for this equilibrium point are 



(Ai,A2)= ( 3(2-7) + V6a,^/- [V6- 



3. Ps ■■ 



and so /C is a source for fc^ < 6 and a saddle otherwise. 
( k I k^\ 

The eigenvalues for this equilibrium point are 

1 



-(fc^-2). 



(Ai,A2 



[6 



[37 



fca] 



(22) 



(23) 



This point exists only for k^ < 6 (when k^ = 6, Ps merges with the equilibrium point K, ). Here, Ps 
is a sink for a < (37 — k'^)/k and a saddle otherwise. 
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4. N : 



37 /37(2-7) + 2aA 
2 A' V 2A2 



kA - 37 



37A: - 2A 



A2 ' ^ A ' 

where A = k + a > 0. Note that this solution is physical (i.e., > 0) either for k^ > 37 or for fc^ < 37 
and 

a> {3"f - k^)/k. (24) 

These solutions were discussed in for a < k and are related to similar power-law solutions discussed 
in This model inflates if 

a>(^7-l)fc. (25) 

(since only a < k was considered in [ p6[ the solutions therein were not inflationary). For fc^ < 2, if 
condition ( p^ is satisfied then (25|) is automatically satisfied and so these models inflate for k^ < 2. 
For 2 < A:^ < 37, if condition (|25) is satisfied then ( p^ ) is automatically satisfied and therefore models 
can infiate for fc^ > 2 if a > (I7 — l)fc. For k^ > 3j there is no constraint on a for the point to 



exist and therefore whether this models inflates is solely determined by (|25|). The eigenvalues for this 
equilibrium point are 



-3 [(2 - 7)fc + 2a] yg [(2 - j)k + 2a]^ ~ 24 [37(2 - 7) + 2aA] [fcA - 37] 

± . . , (26) 



4A 



4A 



and so N is always a sink when it exists. Note that the scalar field acts as a perfect fiuid with an 
equation of state parameter given by 



70 



5. N2 



3(2-7)^ 



7 

2a2 



(27) 



3(2 -72)' 



I/O a 
''=2(^^-2)+ (2^ 



> 0. 



This equilibrium exists for a < W |(2 — 7) and is a saddle, as determined from its eigenvalues: 



(Ai,A2 



:[2-7] 



2a2 



3(2-7)2 



3 a(k + a) 
-7 H — 

2^ (2-7) 



(28) 



Table |^ lists the equilibrium points and their stability for the ranges of k and a. As is evident, the 
presence of the interaction term can substantially change the dynamics of these models. We note that all 
equilibrium points correspond to self-similar models |U0|. 



4 Discussion 
4.1 Inflation 

We first note that we can obtain inflationary solutions when fc^ > 2, unlike the case in which there is no 
interaction term. Moreover, we see from Table 2 that these inflationary solutions, corresponding to the 
equilibrium point N and which occur for a > ^(37 — 2)/c, are sinks (attractors). This result complements 
the results of who looked for inflationary solutions for steeper potentials. When fc^ < 2 and < a < 

(37 — fc2)/2 the power-law inflationary solution corresponding to the equilibrium point Pg is again a global 
attracting solution. 

Of particular interest is the case when k^ < 2 and a > (37 — k'^)/k > 0, when P5 is no longer a sink. 
Therefore, trajectories approach this equilibrium point (and the models inflate for a deflnite but arbitrarily 
large period of time) and eventually asymptote towards the new inflating model corresponding to N. For 
a w {3j — k'^)/k, the eigenvalues for N ( p^ ) are real and negative and so the attracting solution is represented 
by an attracting node. However, for a > (37 — k'^)/2 this equilibrium point is a spiral node; i.e., trajectories 
exhibit a decaying oscillatory behaviour as they asymptote towards N. (For example, for large a equation 
( p6| ) becomes A± w — | ± \/—3ak, leading to complex eigenvalues.) 
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< fc^ < 2 


2 < < 37 


37 < fe"" < 6 


k '^ > 6 




a < Ci 1 a > Ci 


a < Ci 1 Ci < a < C2 1 a > C2 


a < C2 1 a > C2 


a < C2 1 a > C2 


K+ 

(NI) 


R for 2a^ < 3(2 - 7)^ 
s for 2a^ > 3(2 - 7)^ 


(NI) 


R 


s 


Ps 


A 

(I) 


s 

(I) 


A 

(NI) 


S 

(NI) 


DNE 


N 


DNE 


A 

(I) 


DNE 


A 

(NI) 


A 

(I) 


A 

(NI) 


A 

(I) 


A 

(NI) 


A 

(I) 


N2 

(NI) 


s for 2a? < 3(2 - 7)^ 
DNE for 2a^ > 3(2 - 7)^ 



Table 2: The equilibrium points for the model with S = —a(j)^ and their stability for various values of k and 
a. Note that Ci = (87 — k'^)/k and C2 = (87 — 2)k/2. The symbol "R" denotes when the equilibrium point 
is a source (repellor), "s" for when it is a saddle, "A" for when it is a sink (attractor), and "DNE" when 
it does not exist within the particular parameter space. The label "(NI)" denotes non-inflationary models 
whereas "(I)" represents inflationary models. 



An Example 

To illustrate this oscillatory nature, an explicit example is chosen with fc = 1 and 7 = 4/3 (radiation). 
The equilibrium points and their respective eigenvalues are: 



/C±: (Ai,A2)= ±V6a + 2,W- , (29a) 




Ps: (Ai,A2)= U-3,-- , (29b) 



-(3a +1) , V49 + 26a + 330^ - 12a3 

TV: {x,y) = \± = ± ^TTn ' ^^'^ 

2(a + 1) 2(a + 1) 

where a > 3 in order for TV to exist and be a sink and for Ps to be a saddle. Note that for a > 3, N2 
does not exist, /C+ is a saddle and JC~ is a source. Numerical analysis shows that is a spiral source for 
a > 3.65. Figure 1 depicts this phase space for a typical value of a in this range (for illustrative purposes the 
value a = 8 is taken), and the attracting region therein is magnified in figure 2. These figures are typical for 
other values of 7 (this comment is important since we note that in the context of conformally transformed 
scalar-tensor theories, strictly speaking 5 = for 7 = 4/3). 



4.2 Late-time behaviour 

The late-time behaviour of these models, both inflationary and non-inflationary, may also be of cosmological 
interest. For 2 < A:^ < 87 & a < (87 — /c^)//c the late-time attracting equilibrium point is P5 which represents 
a power-law non- inflationary model. Due to recent observations of accelerated expansion models that 

are presently inflating are also of interest. For k^ < 2 Sz a < (87 — k'^)/k the late-time attracting equilibrium 
point is Ps which represents a power-law inflationary model. However, in both of these cases the matter 
contribution is negligible and so these models are not of physical interest. For fc^ < 2 & a > (87 — fc^)/fc 
the late-time attracting equilibrium point is iV, which represents a power-law inflationary model in which 
both the matter and the scalar field are non-negligible and their energy densities are proportional to one 
another. These models are potentially of great significance, and they have been discussed recently within the 
scalar-tensor theory context (see below). Finally, for fc^ > 2 & a > (87 — fc^)/fc the only late-time attracting 
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y 




Figure 1: Phase diagram of the system ^1^ ) when S — —aij)^ for the choice of parameters fc = 1, 7 = 4/3 and 
a = 8. In this figure, the black dot represents the source (i.e., the point ), the large grey dot represents 
the sink (i.e., the point N) and small black dots represent saddle points. The region above the grey dashed 
line represents the inflationary portion of the phase space. Arrows on the trajectory indicate the direction of 
time. 




Figure 2: A magnification of the attracting region of the phase space depicted in figure 1. See also caption 
to figure 1. 
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equilibrium point is N; when a < (87 — 2)k/2 these models represent non-inflating models whereas the 
corresponding models inflate for a > (87 — 2)k/2. 

In the absence of an interaction term, matter scaling solutions are represented by equilibrium points of the 
corresponding dynamical system. We have shown that for simple interaction terms found in the literature, 
these matter scaling solutions cannot be represented by equilibrium points. However, new equilibrium points 
arise which represent solutions in which the energy densities of the matter and scalar field remain a fixed 
proportion to one another and obey 7^ < 7; these solutions are analogues of the matter scaling solutions in 
which 70 = 7 

In | |4l[ | a large class of non-minimally coupled scalar field models with a perfect fluid matter component 
were investigated. These models contain scalar-tensor theory models and, in particular, Brans-Dicke theory 
models with a power-law potential. On performing a conformal rescaling of the metric, the governing 
equations of these models reduce to the equations for a scalar field in general relativity with an exponential 
potential and an extra coupling to the ordinary matter, and are equivalent mathematically to the equations 
studies here (as was noted in [||). Amendola |^ performed a phase space analysis of these models and 
obtained similar results to those obtained here. (Since Amendola assumed < 7 < 2 he obtained a wider 
range of possible behaviours; however, these additional results are of lesser interest in the context of our 
work. In particular, values of 7 > 4/3 lead to negative values for the coupling constant a.) 

The main aim in was to express the solutions back in the original 'Jordan' frame and study the 
cosmological consequences of the underlying scalar-tensor theory models. In particular, 'decaying cosmolog- 
ical constant' solutions were considered which are inflationary and such that the scalar fleld component is 
asymptotically non-negligible. Models consistent with the observations of accelerated expansion and in 
which a physically acceptable fraction of the energy-density is in the scalar field, were found to be severely 
constrained by the upper limits on the variability of the gravitational constant [p3 and by nucleosynthesis 



observations. In further work Amendola |43| considered a quintessential scalar field coupled to matter with 
an additional radiation matter component, and studied the effect of density perturbations on the cosmic 
microwave background in these so-called 'coupled quintessence' models. 

4.3 Early-time behaviour 

From table ^ we can see that the only early-time attractors are /C^ for certain values of k and a, which 
correspond to massless scalar field models (which are analogues of Jacob's vacuum solutions [^). However, 

in table ^ we see that for k^ > 6 and a > ^^(2 — 7) there are no equilibrium points which represent sources 
and the trajectories consequently asymptote into the past towards a heteroclinic cycle. In this cycle, orbits 
quickly shadow the invariant set x'^ + y"^ — 1 {fl — 0), spend a period of time near the saddle IC~ , quickly 
shadow the line y = (V = 0, 0), and then spend an period of time near IC~^ after which the cycle is 
repeated. During each cycle the orbits pass through the infiationary portion of phase space. We stress here 
that this motion is not periodic; on each successive cycle orbits will spend a longer time near the saddles 
JC^ . This past- asymptotic qualitative behaviour, which is depicted in figure 3, is similar to that found in 
within the context of string cosmology; this is not surprising due to the conformal relationship between 
these string models and the models under investigation here [1). 

4.4 The case 7 = 2 

When 7 = 2, corresponding to a stiff perfect fluid, the models is equivalent to a model containing a scalar 
field with an exponential potential and a second interacting massless scalar field. We also note that 7 = 2 is 
a bifurcation value. The equilibrium points and their eigenvalues in this case are 

1. /C+ : (Ai,A2) = |^-\^a,y|[x/6 + fc 
Since a > this equilibrium point is a saddle. 

2. JC- : (Ai,A2) = [\/6a,7|[V6-A: 
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K. is a source for /c^ < 6 and a saddle otherwise. 
3. Ps : (Ai, As) - ( [6 -k'],-[Q~k' - ka] 



Ps is a source for a < (6 — k )/k and a saddle otherwise. 



4. TV: 




n = 



kA-6 



2 [2k 



A2 ' ^ A ' 

where A = fc + a > 0. Note that this solution is physical (i.e., SI > 0) either for fc^ > 6 or for fc^ < 6 
and a > (6 — fc^)/fc. These models inflate for fc^ < 2 as well as for fc^ > 2 and a > 2k. The eigenvalues 
for this equilibrium point are 



A± 



3a ^9a2 - 12aA [fcA - 6] 



± 



(30) 



2A 2A 

and so N is always a sink when it exists and is a spiral sink for a > [Sk]^^ [27 — 8fc^ + y/729 — A8k'^\ 
or a < [8k]-^ [27 - 8fc2 - V729 - 48fc2J (and is a spiral sink for all fc^ > 6). Note that the scalar field 
acts as a perfect fluid with an equation of state parameter given by 

2 



70 



1 + iaA 



<7- 



(31) 



Note that the equilibrium point N2 does not exist for 7 = 2. Table |^ summarizes the stability analysis 
for 7 = 2. The qualitative behaviour is similar to that of the case 7^2, except that there is no region 
corresponding to 37 < fc^ < 6 when 7 = 2. Note again that there exists a heteroclinic cycle at early times 
for k > \/6 (see Fig. 3); unlike the general case 7 7^ 2, this heteroclinic cycle exists for all a > 0. 





< fe^ < 2 


2 < fc^ < 6 


fe^ > 6 




a < Ci 1 a > Ci 


a < Ci 1 Ci < a <2k \ a > 2k 


a < 2k 1 a > 2k 


(NI) 


s 


ic- 

(NI) 


R 


s 


Ps 


A 

(I) 


S 

(I) 


A 

(NI) 


S 

(NI) 


DNE 


N 


DNE 


A 

(I) 


DNE 


A 

(NI) 


A 

(I) 


A 

(NI) 


A 

(I) 



Table 3: The equilibrium points for the model with 5 = —a(j)^ and their stability for various values of k and 
a with 7 = 2. Note that Ci = (6 — k'^)/k. See caption for table ^for notation. 



5 Interaction term of the form 6 = afiH 

This section provides a second example to demonstrate that other types of interaction terms can also lead to 
similar behaviour; i.e., that there will be a range of parameters for which the inflationary models which drive 
the matter fields to zero are not late-time attractors and for which the trajectories exhibit an oscillatory 
behaviour as they asymptote toward the late-time attracting solution. Specifically, the interaction term 
5 = a/iiJ is chosen, where a > 0. 

With this choice, equations dill ) become 



-3a;(l - x^) 



y 



y^fcx -I- 2a;2 + 7(1 



(32a) 
(32b) 
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y 
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For physical reasons we are not interested in early-time behaviour, and hence the line a; = will not be 
considered. Consequently, a full phase-plane analysis is not possible using these variables. However, it is still 
possible to determine the equilibrium points with x ^ Q for the system and determine their local stability. 
There are four equilibrium points for this system for x ^ Q: 

1. /C+: Eigenvalues for this equilibrium point are 

(Ai,A2)- |^3[2-7]-|-a,y|x/6 + 3j (33) 

This equilibrium point is a source. 

2. /C^: Eigenvalues for this equilibrium point are 

(Ai,A2)= (^3[2-7]+a,-y| + 3y (34) 

K.~ is a source for fc^ < 6 and a saddle otherwise. 

3. Ps- Eigenvalues for this equilibrium point are 

(Ai,A2)= (^-\[Q~e],~[S-k^-i-k^]^, (35) 

where 5 = 37 — a. This point exists only for < 6, and is a source for k^ < 37 and a < 37 — /c^ 
(saddle otherwise). 



Note that this solution is only physical (il > 0) for [37 — k^] < a < 3j and k^ < 37, and represents 
an inflationary model for a > (37 — 2) (consequently, this model will always inflate for k^ < 2 and can 
inflate for 2 < fc^ < 37). The eigenvalues for this equilibrium point are 



A± 



[2a(2S-P) + (2-7)S^] 
47S 



/[2a(2S - fc2) + (2 - 7)S2]^ F „ 8jE^{k'^ ^ E) [2aP + (2 _ 
^ 47AA: ' 

(36) 

and so Ni is a sink for a < 37 — ^fc^. Note that the scalar field acts as a perfect fluid with an equation 
of state parameter given by 

= TT^ ^ ^' ^^^^ 

Table ^ lists the equilibrium points and their stability for the ranges of k and a. Again, for the range 
k^ < 2 and (37 — ■^k'^) < a < (87 — |fc^), the power-law model Pg is no longer a sink and Ni is a source. 
Therefore, solutions approach the equilibrium point which is represented by Pg (thereby inflating) for an 
indefinite period of time, but eventually evolve away. It can be shown numerically that within this range 
for k and a, the equilibrium point A'^i is a spiral node (for instance, for k = 1 and 7 — 4/3, A^i is a sink for 
3 < a < 3.5 and a spiral sink for 3.24 ^ a < 3.5); therefore the scalar field exhibits oscillatory motion as the 
solutions asymptote toward N. Care must be taken in interpreting the analysis and obtaining global results 
since this system is not well defined for x — {(j) ^ 0) for this particular example. 
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< A:^ < 2 


2 < /c^ < 37 


37 < A;^ < 6 




a < C3 


a > C3 


a < C3 


a > C3 








a < C4 


a > C4 




a <Ci 


a > C4 












R 

(NI) 










A 


s 


A 




s 






(I) 


(I) 


(NI) 




(NI) 




TV 




A 


s 




A 


s 






DNE 


(I) 


(I) 


DNE 


(I) for 

> 3t - 2 


(I) tor 
a > 37 - 2 


DNE 



Table 4: T/ie equilibrium points for the model with 5 = afiH and their stability for various values of k and 
a. Note that C3 = 37 — and C4 = 37 — ■^k^. The symbol "R" denotes when the equilibrium point is a 
source (repellor), "s" for when it is a saddle, "A" for when it is a sink (attractor) , and "DNE" when it does 
not exist within the particular parameter space. The label "(NI)" denotes non-inflationary models whereas 
"(I)" represents inflationary models. For fc^ > 6 the only equilibrium points to exist are the points JC^ ; JC~^ 
is a source and 1C~ is a saddle (i.e., there is no sink). Recall that the system is not well defined for x = 
and therefore a global analysis is not possible. 

6 Conclusions 

Without an interaction term, it is known that for k^ < 2 the global late-time attractor for the system ( [III ) 
is a power-law inflationary model in which the matter is driven to zero ]3^ ]. The purpose of this paper was 
to show that this behaviour could be altered qualitatively with the introduction of an interaction term. In 
particular, for models with an interaction term of the form 5 — —acpfi there are values in the parameter 
space for which the equilibrium point Pg, corresponding to this particular power-law inflationary model, 
becomes a saddle, and so while the models may spend an arbitrarily long period of time inflating with 
^ 0, they eventually evolve away from this solution. The late-time attractors within the same parameter 
space, corresponding to the new attracting equilibrium point A'^i, are also inflationary but with the matter 
field's energy density remaining a fixed fraction of the scalar field's energy-density and with 7^ < 7. These 
are analogues of the matter scaling solutions in which 7^ = 7. 

For an appropriate parameter range, the equilibrium point A^i is an attracting focus, and hence as 
solutions approach this late-time attractor the scalar field oscillates. Although the late-time behaviour (cor- 
responding to Ni ) is still infiationary, the oscillatory behaviour provides a possible mechanism for infiation 
to stop and for conventional reheating to ensue (indeed this is similar to the mechanism for reheating in 
scalar field models with a potential containing a global minimum p^-|2^ ) . To study reheating properly more 
complicated physics in which the oscillating scalar field is coupled to both fermionic and bosonic fields needs 
to be included. This contrasts with the situation for exponential models for k^ < 2 with no interaction term 
which have no graceful exit from infiation and in which there is no conventional reheating mechanism. 

Therefore, we have shown that there are general relativistic scalar field models with an exponential 
potential which evolve towards an inflationary state in which the matter is not driven to zero and which 
exhibit late-time oscillatory behaviour; these models may constitute a first step towards a more realistic 
model. There is the question of how physical these models are, since they correspond to relatively large 
values of a. In the context that the interaction term represents energy transfer, for physical reasons it might 
be expected that a must be small; i.e., a < 0.1 ||3^ (see also ^]). On the other hand, in the context 
of scalar-tensor theories a is of order unity and can certainly attain values large enough to produce the 
behaviour described above (see Eq. (||); this is also the situation in the context of string theories. 

It is also of interest to study the cosmological consequences of the 'decaying cosmological constant' 
or 'quintessential' cosmological models, since they may be consistent with the observations of accelerated 
expansion 0, D and may lead to a physically interesting current residual scalar field energy-density. These 
issues have recently been addressed by Amendola [^,|4^ in the context of the conformally related scalar- 
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tensor theories of gravity. 
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